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Hamilton's equations with noise and friction possess a hidden supersymmetry, 
valid for time-independent as well as periodically time-dependent systems. It is 
used to derive topological properties of critical points and periodic trajectories in an 
elementary way. From a more practical point of view, the formalism provides new 
tools to study the reaction paths in systems with separated time scales. A 'reduced 
current' which contains the relevant part of the phase space probability current is 
introduced, together with strategies for its computation. 



I. INTRODUCTION 

Morse theory describes the relationship between the critical points of a smooth real 
('Morse') function and the topology of the manifold on which it is defined [if. It is a 
major tool in mathematics, but also finds natural applications in physics, for instance in the 
classification of periodic orbits in classical mechanics. 

An elementary and elegant derivation of Morse theory was obtained years ago P] through 
the use of supersymmetric quantum mechanics (SUSY-QM), in which the potential is related 
to the Morse function. It is based on the fact that the semi-classical lowest eigenstates of 
the SUSY-QM Hamiltonian are concentrated on the saddle-points of the function - those 
having k fermions on saddles with k unstable directions. The supersymmetry operators that 
map eigenstates with k fermions into eigenstates with k ± 1 fermions then induce relations 
between the corresponding saddles. 
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Soon after this development, it became clear that the SUSY-QM Hamiltonian - or, more 
precisely, its zero-fermion restriction - is related by a change of basis to the (Fokker-Planck) 
equation describing the evolution of probability associated with a Langevin process 0, % 
1^, 0], the semi-classical limit now becoming the low temperature limit. In this basis the 
lowest eigenstates with k fermions are concentrated, rather than on saddle points with k 
unstable directions, on the unstable manifolds emanating from them. For example, one- 
fermion eigenvectors with small eigenvalues are peaked on the gradient paths joining two 
minima via a saddle point, and so represent the reaction current jjjj. 

A low-temperature Langevin process is a practical method to locate minima of the po- 
tential: the dynamics consist of gradient descents, and the small noise allows to escape local 
minima - a form of what is known as 'simulated annealing'. In the language of SUSY-QM, 
such a process corresponds to the evolution in zero-fermion subspace. One is then led to 
ask whether a dynamics associated with the one-fermion subspace can give a 'simulated 
annealing' scheme that will converge to reaction paths, just as the ordinary one does for 
metastable states. Indeed, this is so [3,0], and has been proposed as a basis for numerical 
algorithms. 

In practical applications it is sometimes necessary to extend the overdamped Langevin 
treatment to a case where inertia plays a role. This is possible because pure Hamilton's 
equations have themselves a supersymmetry, whose consequences have been explored ex- 
tensively by Gozzi, Niemi and their co-workers . There is however a 
problem that makes the formalism less transparent than in the Langevin case: in classical 
mechanics the (Liouville) equations for the evolution of probability in phase-space do not 
contain second derivatives: the spectrum of the evolution operator is then continuous and 
the space of wavefunctions awkward. A natural cure to this problem, both from the ana- 
lytical and the practical point of view, is the introduction of noise and friction so that we 
study a more general process: 

I q%=p% 

V(q) is a potential energy, rji are Gaussian white noises modelling the interaction with a 
thermal bath at temperature T, while 7 is the coupling to the bath (physically a friction 
coefficient). The dynamics (JTJ) can be expressed as a probability density evolving according 
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to the Kramers equation 



dPjq, p,t) 
dt 




) 



d_ 



Pi P(q,p,t) 



H K P{q,p,t). (2) 



It turns out that one can uncover a hidden supersymmetry associated with this equation, 
just as there is SUSY-QM associated with the Fokker-Planck equation without inertia. As 
we shall see, this leads to a non-Hermitian supersymmetric quantum mechanics, whose zero- 
fermion restriction is the Kramers equation. The higher fermion-number subspaces contain, 
again just as in SUSY-QM, the information on the fixed points and their stable and unstable 
manifolds, where stability is now defined with respect to Hamiltonian dynamics perturbed 
by friction. Perhaps more surprising is the fact that one can generalize these results to the 
case in which the Hamiltonian depends periodically on time - one then has a supersymmetric 
structure in the Floquet representation. This supersymmetry has a series of consequences 
which we shall also explore. 

The organization of this paper is as follows: in Section II we construct the extension 
of the Kramers operator and we discuss the consequences its supersymmetry has on the 
organization of the eigenvectors and eigenvalues. In Section III we study the low-temperature 
('semi-classical') limit using standard path-integral methods, and show that it is dominated 
by periodic trajectories - in particular fixed points - of the noiseless dynamics. In Section 
IV we use these results to show how one can rederive Morse theory results using time- 
independent Hamiltonians. As an alternative, a WKB treatment of the time-independent 
case with conservative forces is possible (Section V): it does not rely on path integrals and 
is very close to the treatment of SUSY-QM, with the only complication of non-Hermiticity. 
In section VI we present the supersymmetries associated with a periodically time-dependent 
Hamiltonian, and give a first few applications. 

Just as the zero fermion subspace corresponds to Kramers' equation, and this in turn 
to a process following Hamilton's equations plus noise and friction, one may ask to what 
stochastic processes does the k > fermion subspaces correspond. In section VII we devise 
such processes, and use them to give a constructive (and quite non-rigorous) derivation of 
the low-temperature wavefunctions yielding the Morse complex. As mentioned above, part 
of the motivation for this work is the construction of algorithms to find saddle points and re- 
action paths between metastable states, a very important problem in Physical Chemistry. In 
Section VIII we show the relation between one-fermion wavefunctions and reaction currents. 
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Interestingly enough, the formalism strongly suggests that rather than studying the currents 
themselves, a modified 'reduced current' should be used, which contains only the part of the 
current that is effective in making transitions, equilibrium circulations within states being 
subtracted. Finally, in the Conclusion we outline the several possible continuations of this 
work. 



II. SUSY OF KRAMERS DYNAMICS 



A. The Hamiltonian case 



The Hamiltonian dynamics 

. an . m ... 

where 7i is the Hamiltonian function, induces an evolution for probabilities P(q, p) in phase- 
space given by 

— --HP with H - VV— — - — — ^1 U) 

dt V dq { dpi dpi dqi ) ' 

One can uncover a group of symmetries of by extending the space with AN fermion 



operators (oj, aj, bi, b\), and writing |9 
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„ s _^(dnd dnd\ d 2 n hj d 2 n t d 2 n t t 

H ~ trl%% d Pl d qi ) + d Pt d Pj a ^ + d Pj d qi {b ^ a i aih {b) 

which reduces to the original in the zero-fermion subspace. has a large group of 
symmetries, generated by the operator whose action is to multiply by 7i (and other constants 
of motion, if present), and by 

N N N 



i=l i=l i=l 

" / a a 



(6) 
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The supersymmetric charges are nilpotent: 

Ql = Ql = Ql = Ql = o, (7) 

and 

K = (Qi + Q 2 + Q 3 ) 2 - (8) 

As already mentioned in the Introduction, the question of the spectra of operators and 
the underlying Hilbert space is rather tricky, because all operators have at most first order 
derivatives. Here we shall work with the Kramers equation, for which all relevant operators 
have discrete spectra, the Hilbert space is tractable, but the symmetry group is considerably 
smaller. 



B. Extended operator and symmetries. 

Let us go back to the original Kramers dynamics with inertia and dissipation (0), which 
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we shall write in a more general form 

OH 

dqt dpi 

Here 7i is a general Hamiltonian function of {q, p} coordinates, equation (0) corresponds 
to H, of the particular form 7i = p 2 /2 + V(q). We shall hence assume that 7i is a smooth 
function, either defined over a bounded phase-space {q, p} or growing fast enough at infinity. 
The probability density in the phase space evolves as: 
dP(q,p,t) 



Of 



-H K P(q, p,t). 

with h =- flfr- ^Vff-----1 (10) 

7 i^x V ®Pi) ~t v ^<li dpi dpi dqi J 

One can easily see that the Gibbs density (e~^) is the stationary state of the process, the 
eigenstate of with zero eigenvalue. 

In order to construct a fermionic extension we start from the observation that the previous 
symmetry charge Q\ is independent of the dynamics considered, having only a geometrical 
meaning (it is related to the so-called exterior derivative). We then propose 

N f f) f) \ 

= (^+^)> (") 



i=l 
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as one of the supersymmetric charges. By inspection, one can see that 

N 



id — 1 

(12) 

has Q as a symmetry. By analogy with the Hamiltonian case one may ask if there is a first 
order differential operator Q satisfying: 

g 2 = Q 2 = 0, T[Q,Q] + =T(Q + Q) 2 = H. (13) 



This is so with: 

N 



* .v^Lt/^ 10H\ f fd ldH\ L+ / d ldH 



(14) 
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i=l 

Thus, we now have a supersymmetric extension of Hk and the corresponding charges 
The large algebra of symmetries (JBJ) of the classical formalism has become much smaller; 
just as in the case of SUSY-QM, only Q, Q and F are symmetries. 

The system (J9j) can be written in terms of the variables ...,X2n) — (QjP) m a more 
compact way as: 

* = -tv.g - d„ - ^ o y = (- ;■;) B „, = 7 (; :). (15) 

When Z)jj = 0, (|T5|) is the usual simplectic formulation of Hamilton's equation, the second 
term of the r.h.s represents the interaction with the bath. With this notation, the generalized 
Fokker-Planck Hamiltonian (|12|) becomes: 

= H K + A ik c\c k , 
where (ci, C2n) = (ai, «tv, &i, &jv), and we have defined: 

A * = (c « + ^aSi/ (17) 

The charges are: 

27V Q 
O 



dxi 

i=i 



(18) 
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The notation in this section has been specialized to the case in which all the phase space 
velocities derive from a global function 7i. Let us keep in mind, however, that when this is 
not the case, and the dynamics is given by: 

±i = -SlijHj - DijHj - J —r)i, (19) 

with 1^ = but Jijix) not globally a gradient, there still is a supersymmetry: 

2N . 2N , r, \ 

- + H J )c\ 

(20) 



i=l 1 i=l v J 



H = ~llx \ Dij ( T ^" + H *) + QljH ^ + Ai 'J^ Ck With Aik = ( Di 3 + ^ 



3. 

C. The spectrum of H 

We shall first examine the case of a time- independent Hamiltonian Tl, including cases like 
(|20p where there is no global potential. Section fVTl will be devoted to the time-dependent 
case while sections IIVI and |V] will be specifically dedicated to the conservative case with a 
global Tl. 

H is not hermitian, and cannot be taken to a Hermitian form via a similarity transfor- 
mation. It acts on functions of the form 

IV') = £ ^.....^....jnCq, P)bl ■ ■ ■ b\J n . . . a]j-), (21) 

where |— ) is the Fermion vacuum. We have to distinguish right and left eigenvectors. They 
are defined as: 

H\lf) = Xili/f) and Ht\tf) = %\tf). (22) 

In what follows we shall suppose that H is diagonalizable i.e. one can find a bi-orthonormal 
eigenbasis 2s| : 

(^f} = ^ lV (23) 

The purpose of the two following subsections is to show that the organization of the 
spectrum is as in figure [TJ 
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FIG. 1: Generic representation of the spectrum of H. The eigenstates are divided in two kinds: 
pairs of eigenstates made via Q and unpaired eigenstates, which does not belong to any pair by Q. 
The latter are only present in the A = eigenspace. 

1. Eigenvectors with non-zero eigenvalues 

Eigenvectors with A ^ are not annihilated by both Q and Q. As Q and Q commute 
with H, they map eigenvectors into eigenvectors with the same eigenvalue. This yields for 
the A spectrum a structure of degenerate pairs connected by Q and Q. One can indeed 
construct a basis (|0f ), |xf)) such that: 

Wf> = lx?> and g|xf> = ^l0f>. (24) 

To see this, first notice that a general eigenstate of eigenvalue A ^ is the sum of 
two eigenstates \<fi R ) and \x R ) annihilated by QQ and QQ respectively. Indeed, denoting 
\X R ) = QQ\^ R ) and \(f) R ) = QQ\4) R ) we have: 

\^ R ) = \h\^ R ) = j (\X R ) + • (25) 

Using the fact that QQ and QQ both commute with H, \x R ) and \<f> R ) are also eigenvectors of 
H with eigenvalue A, and are annihilated respectively by QQ and QQ. One thus constructs 
a basis of the whole eigenspace A ^ as the union of bases of eigenvectors annihilated by 
QQ on one hand, and QQ on the other hand. We do this as follows: denote |xf ) a basis of 
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eigenvectors annihilated by QQ, and define = j-Q\x R )- Because Q\4>f) = -^H\xf) = 
\xf)i the \4>f ) so defined are independent, since one can map them back into an independent 
set with Q. 

Furthermore, the \(j) R ) generate all the eigenvectors with A ^ annihilated by QQ. 
Indeed, take one such \ip R ). One clearly has \ip R ) = \QQ[^p R )- As Q\ip R ) is annihilated 
by QQ, it can be developed as Q\ip R ) = X^^lxf); an d hence \ip R ) = ? JZ f a iQ\xf) = 

The family (|0f ), |xf)) is thus a basis of the whole Aj 7^ spectrum, which satisfy the 
pairing property (j2U). As to the left eigenvectors, the structure is the same. Constructing a 
basis with 

(x-\xf) = ^ 

(26) 

(<f>f\xf) = (xM) = °> 

the left eigenvectors are paired according to: 

QM) = \ti) QW> = ^lxf>- (27) 



2. Zero eigenvalues and topology 

The scenario in the A = eigenspace is a little more complex. Let us argue that one can 
in principle build a basis of right eigenvectors (we drop the index R within this subsection) 
composed of (see figure HJ): 

i) pairs (|0f +1 ), \xD) such that \ X \) = Q\$ +1 ) ± 

ii) unpaired eigenstates \p\) such that Q\p\) = and V|^>) \p\) 7^ Q\ip), 

where k denotes the number of Fermions. Note that such a basis is matched by the corre- 
sponding one for the left eigenvectors. 

We shall now construct a basis satisfying ()28|). Let us first look at the fermion sector. 
All the eigenvectors are annihilated by the aj and consequently by Q. Some of them are the 
image by Q of other eigenvectors, and some are not. Let us note \x®) a basis of the former, 
and the 1-fermion eigenvectors which generate them: \Xi) = Q\4>})- The Ix®) may not 
be a basis of the zero fermion sector and one can complete them with eigenvectors |p°). By 
definition, the \p° i ) are annihilated by Q but are not the image of any other eigenvectors 
by Q. At this point, \p°i)} constitutes a basis of the fermion sector satisfying (|2~%|) . 
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One can then turn to the 1 fermion sector. The part of this sector which is annihilated 
by Q can be organized as the fermion number, that is there exist a basis {Ix*), I Pi)} 
which generates this part, and eigenvectors \4>f) such that \xj) = Q\<Pl)- The family \(f)\) 
introduced above completes the {|Xi)> |0i)} i n a basis of the whole 1 fermion sector. This 
construction can be followed for every fermion sector and allows us to construct inductively 
a basis {\Xi), \4>i), \Pi)} which satisfy (|25j). 

In the usual SUSY-QM, only unpaired states exist in the zero eigenvalue subspace. Indeed, 
we shall show in Section IS that this is also the case for a Kramers problem with time- 
independent forces that derive from a global potential. Note however, that states paired by 
Q with A = can exist: a system living on a ring encircling a magnetic flux has no global 
potential and there is a pair by Q. 

The eigenstates |pf) are related to the topology of the phase-space because is the 
exterior derivative acting on (left) states with k fermions (the differential /c-forms). They 
span a space whose dimension is the k th Betti number B^, and is isomorphic to the so-called 
k th de Rham cohomology group associated with the phase space. 

The basis of the A = and A ^ eigenspaces form a global basis: {|^i)} = 
{|xf), |0f), |pf)}- Given its structure, it is tempting to denote the states generated by 
{|pf )} an d {|0f) ? Ixf)} as "unpaired" and "paired", respectively, even if one can construct 
states which are not paired by Q without being generated by the \p\). 



III. LOW-TEMPERATURE LIMIT: FIXED POINTS AND PERIODIC ORBITS 

A. Fixed points and periodic orbits 

To study the low-lying eigenvectors in the limit of small T, we shall compute the trace of 
the evolution operator for different values of t and use the result to reconstruct the spectrum. 
A quick method to do this is to write a path integral and use saddle point evaluation. This 
is a standard exercise which we report in Appendix El When the temperature goes to zero, 
one finds that the path integral is dominated by the noiseless periodic trajectories which 
satisfy: 

on m 

oqi dpi . , 

, c m (29) 
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To obtain the next order, one develops the coordinates as a small perturbation x\ = {q'^p'j) 



it, Pi) 



around each noiseless orbits x\ ■ 

Xi = x\ + VTx'i. 

The contribution of each orbit x c can be seen as the trace of Te~ fo # c (*') d *' with: 



(30) 



d 2 n 



dpidpj 



d 2 n 



dpidqj 











I 7-^ 1" IP* 

dp'i \ dpi 3 



d 2 n 



dpidpj 



Q ' ,P 



+iq 3 

+ 7 



d 2 n 



dpidqj 
d 2 H 



i / 
+ Pi 



d 2 n 



Q ' ,P 



dpidpj 



3 dqidpj 

d 2 n 



+ qj 



, d 2 n 



dpidp 



a% + 



dqidqj 
d 2 H 



+ 



d 2 n 



dqidqj 



b\a.j 



9 tP 

(76] at + b\bj -a)ai) 



(31) 



1 iP 



k \ 7 / 



where A^(t) is defined as in (|17|). but evaluated along the classical periodic orbit (which in 
certain cases will just be a fixed point): 

&H 



AUt) = A k [x c (t)} = (Ai + %) 



' dxkdxj 



(32) 



Here we have adopted a notation involving a global potential 7i, although this is not neces- 
sary: if there is no global H it suffices to write (}3~T]) with A c k - defined as in (J20J). 

Note that to this order fermionic and bosonic parts are decoupled. (13 1J) is nothing but 
the SUSY Hamiltonian corresponding to a diffusion on a time-dependent harmonic potential 



7i c = \ qx-qL x i x p corres P on ding to 



x, 



(33) 



-A%{t)x'j + DiTji, 

where rji is a Gaussian white noise. 

A compact way of expressing the spectral properties of the evolution operator is via the 
generating function: 



T(\,t) = Tr (\ F Te-5 H ^ dt ') = ]T A fc Tr (r e ^ H ^ dt ' 



k fcrm. 



(34) 



where T denotes time order and F is the Fermion number (JBJ). To leading order in the 
temperature T, we have to compute T(X,t) as a sum over the contributions T c (X,t) around 
each periodic orbit: 



T(A,t)~VA fc V Tr(Te-/"W) ee V T%X,t) 

J V / k ferm. J 



noiseless 
orbits c 



(35) 



noiseless 
orbits c 
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Because boson and fermion degrees of freedom are decoupled to this order, each term of 
the sum is a product of a trace over boson and a trace over fermion degrees of freedom 
- both with a (in general time-dependent) harmonic oscillator. Again, this is a standard 
exercise which we do in detail in Appendix El The result is that for every orbit the boson 
degrees of freedom contribute with a factor | det(l — ?7 c (t))| _1 , while the fermion ones with 
det(l + XU c (t)), where U c {t) is a 2N x 2N matrix defined by: 

U c (t') = -A c (t')U c (t') U c (0) = 1. (36) 

We have then, to leading order in T: 

v det ( 1 + A ^)) (37) 

noiseless 1 v w/l 
orbits c 

This formula has two limitations, both reflecting important features of the phase-space 
structure: 

• If an eigenvalue of U c is 1 for some orbits, (|37|l diverges. This can be accidental, e.g. a 
critical point that is undergoing a second order phase transition, or, more importantly, 
a consequence of the fact that the orbit is not isolated but belongs to a continuous 
family, possibly as a result of a symmetry: the problem becomes one of degenerate 
Morse theory 26j . 

• The saddle point evaluation is legitimate to the extent that T — > at fixed t. If we are 
interested in orbits of period going to infinity as T — > 0, we have to bear in mind that 
the action itself will depend on T, and (|3T|) may be invalid. This problem is not specific 
to our treatment : the long time periodic orbits are a usual pitfall in semi-classical 
quantization. 

B. Spectrum of low real eigenvalues 

Let us analyze first the contribution to (|3T|) of an isolated orbit whose period equals the 
period r of the Hamiltonian: 

HA,,r)~nW. ,38) 

i=l 1 V V I 
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where (%,..., w 2 at) are the eigenvalues of U c (t) (for simplicity, we omit the supra-index c 
labelling the noiseless trajectory). One can read in the factor (1 + Auj) the contribution of 
the Fermion vacuum whose eigenvalue is 1 plus A times the contribution of the one-fermion 
state, whose eigenvalue is U{. The bosonic counterpart can be read in -, : we shall expand 
\i2 u r}\ as a geometric serie from which we shall recognize the spectrum. 
First consider the case of a real U{. If \uA < 1, then 
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I 1 

^ = ^^ = Y<\ (39) 
If \u,i\ > 1, one can rewrite |1 — Ui\ as |w.t|(l — ^-). The same development gives: 

-— — = sign(n,) "7* V*>1. (40) 
' U{ ' fc=i 

Next, consider the case of a pair of complex eigenvalues (ui,u*). If \ui\ < 1 we may write: 

£ ( Ul r k (u;r k \ (4i) 



(1-<)(1-M* n ) 
ill « I V * 7 fc=0,fc'=0 



while, if > 1: 

1 1 = E (u t )- nk (u*y nk '. (42) 

fe=l fc'=l 



m! 1 ! II — u* 



All in all, the contribution of the orbit to the spectrum is the tensor product of the following 
sets (see figure EJ): 

• (1, Ui) <g> (1, Ui, uf, ...) for Ui = u*, \ui\ < 1; 

• (!> u i) ® ^> •••) for M i = w *> kl > 

• <g> (l,Mi,M-, ...) <g> (1,<) <g> (l,u*,uf ,...) for Mi ^u*, \Ui\ < 1; 

• (1, Ui) ® (i, ^, ...) ® (1, <) ® (i, ^, ...) for Ui ^ <, |ui| > 1. 

We have supposed here that the number of real eigenvalues U{ < — 1 is even. 

Let us now consider the contribution of an orbit of primitive period pr. Clearly, one can 
start from p different points along the orbit. For each starting point the preceding discussion 
holds, and one gets for the evolution over time pr a spectrum as above, but each level is 
now p-fold degenerate (to this order in T). The path integral tells us that the contribution 
of this orbit to any trace over n cycles is zero if n is not a multiple of p, and is p times the 
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\Ui\ < 1 



|A| 



1 

\ u i\ 



|A| 



1 

\ u i\ 



Wj > 1 



I"/ 
1 



|A| t 



\Ui\ 



Y 

Fermionic 



Y 

Bosonic 



FIG. 2: Structure of the spectrum for \ui\ < 1 or \m\ > 1 at leading order in T. The global 
contribution of the orbit is the tensor product of the bosonic and the fermionic parts. There is a 
gap of order one in modulus between the most stable state and the first excited one when T goes 
to zero. 



contribution of a single starting point otherwise. This can be understood if each multiplet 
yields, for the evolution over a single period of the Hamiltonian, eigenvalues corresponding 
to the p different p th roots of those for p cycles — because the sum of the p different roots 
to the power n is non-zero only if p divides n. 

Let us stress that breakdown of (|37[). for example because one of the \ui\ = 1, signals 
the fact that the spectrum is no longer a superposition of harmonic spectra of frequencies 
of order one, and that the gap between ground and first excited state will go to zero with 
T. In such cases, the contribution of the it, such that \ui\ 7^ 1 is as described here, and the 
degrees of freedom in the directions corresponding to \uj\ = 1 have to be treated with other 
methods (collective coordinates, for example). 

The main point of this section is that each orbit is associated with one and only one 
eigenstate of the evolution operator with unit eigenvalue to this order. This eigenstate has 
k fermions if U c (t) has k eigenvalues with modulus larger than 1: that is to say if its Morse 
index is k. 



15 



IV. TIME-INDEPENDENT HAMILTONIAN: A NON-HERMITIAN SUSY 
QUANTUM MECHANICS AND MORSE THEORY. 

In the case of a time-independent Hamiltonian deriving from a global potential, we recover 
the Morse inequalities for the stationary points of the dynamics, by going to the small 
temperature limit. The geometric structures involved now correspond to the manifolds that 
are stable and unstable with respect to Hamiltonian dynamics plus friction, instead of simple 
gradient descents. 

If 7i is time-independent and conservative, in the presence of non-zero friction, the only 
periodic orbits that matter are fixed points, as one can see using 



dH ^dH. c dH. c ^fdH\ 2 

Because the energy along an orbit has to be periodic, the only possibility when 7 7^ is that 
it is constant and |^ = 0. If the Hamiltonian is of the form Ti = \p 2 + V(q), this implies 
that Pi=0 and ql = constant. More gene rally , this implication has to be verified, but it is 
true in all but very pathological examples |27j]. 

For every fixed point x c , the eigenvalues u\ of U c (t) are 



< = e"^, (44) 

where the are the (in general complex) eigenvalues of Aj 3 -[a; c ] (Cf. equations (}32|) and 
(136)1 ) . The Morse index of such a critical point is the number of eigenvalues such that \u^\ > 1 
or equivalently ReA^ < 0. The result of the preceeding section implies on one hand that all 
the moduli of the eigenvalues of e~ tH are to this order in T smaller or equal than one. On 
the other hand, the number of eigenvalues that are one within the k fermion subspace 
coincides with the number of critical points of index k. 
For large t, we have then: 

gnT(A,t) = Hm£A*Tr (OU = E A " M *' ( 45 ) 
k k 

because eigenvalues of e~ tH with moduli smaller than one are exponentially suppressed as 
in ordinary SUSY-QM. 
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On the other hand, 

27V 

T(X,t) = J^A fc Tr (e""*) 



k fcrm. 



k =° (46) 

27V 27V v ' 

= y, (--'") irrr + e ^ irsi . 

fc=0 fc=0 

where the supraindices denote traces taken over subspaces spanned by "paired" and "un- 
paired" states, as defined in section III C 21 We now introduce the partial traces: 

D / % / I paired states anihilated by Q f Ai-r\ 

K k [t)=li{e j| fcfcrm . (47 ) 

Taking into account the pairing of the spectrum (figure Q}, the fact that the eigenvalues of 
the unpaired eigenstates are zero and that the dimension of the space they generate gives 
the Betti numbers (see Section Hi C2|) . we get: 

27V 27V 



T(A, t) = Y, xkBk + E xk (Mt) + Rk-i(t)) 

27V 27V ^ > 

hm T(A, t) = ^2 \ k B k + J2 Xk iMoo) + Hfc-i(oo)) 



t— >oo 

fc=0 fe=0 



where the R k (oo) are integers: the number of paired eigenstates of H annihilated by Q 
having eigenvalue zero to leading order in T. Putting together (}45|) and (jUJ), we have: 



M fc = B k + i? fc (oo) + fl fc _i(oo). (49) 

The positivity of the Rk(oo) (except for i?_i(oo) = 0) constitute the strong Morse inequali- 
ties: 

v p 
Vp ]T(-l) fc iV fc < ]T(-l) fc M p _fc (50) 

fc=0 fe=0 

It is an easy calculation to show for Tt of the form + V(x) that the index defined 
here as the number of eigenvalues with negative real parts tends to the usual index defined 
as the number of negative eigenvalues of the potential V as 7 — > 0. 

Before concluding this section, let us remark that all we have done here is valid if 7 is 
kept finite as T — > 0. If we wish to consider 7 — > together with T = 0, then orbits that are 
not fixed points contribute. This is most easily seen by considering the action in the path 
integral of Appendix El which can be written: 

1 r J /. on on\ 2 1 r , / on\ 2 7 r J . 2 [, 

5 = — - dt (pi + 7- — = — - / dt Ui - — + / dt q£ - dt — , 



(51) 
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restricted to periodic trajectories such that qi = If the Hamiltonian is time-independent, 
the term ^ J dt q 2 suppresses periodic orbits that are not fixed points. However, if we 
consider 7 — > 0, T — > 0, keeping ^ finite, then non-fixed orbits will contribute, and we may 
get richer inequalities. We shall not follow this strategy here, although it seems a promising 
line of research. 



V. A WKB APPROACH. 

In this section we treat a Hamiltonian of the form: 

%p) = ^| + y(g), (52) 

i 

with a standard WKB treatment that allows to derive Morse theory with a construction 
very close to the one in SUSY-QM. 



A. Conservative forces 



For Langevin processes (without inertia), the Fokker-Plank operator is non-Hermitian but 



can be brought to a Hermitian form by a symmetric real transformation 3, 15], provided 
the forces are conservative. This is true also for the extension with fermions, and the 
transformation yields the Hermitian SUSY-QM. In the Kramers case H cannot be brought 
to a Hermitian form, and part of the difficulty with the proofs comes from there. In this 
section we shall take a closer look at a case for which the discussion simplifies considerably, 
that of Hamiltonians of the form (15211 for which H reads: 



H 



N 

E 

i=l 



d 2 dV d d d 

oPi oqi 



lT dp 2 dqidpi 



N 



d 2 V 



— i 



N 

i=i 



(53) 
(54) 



One can easily verify that H, Q and Q are then related to their adjoint by [28]: 

H ] = RHR~\ RQR- 1 = Q\ RQR' 1 = Q\ 

where R is a real, Hermitian invertible operator defined as (no summation inside the square 
brackets) : 

R = R^ = e^PJ, where P|^(q, p)) = |-0(q, -p)), 
J = JJ[1 + a\bi + b\ai - a\ai - b\bi + ^a\ai + (7 - 2)a\b\aibi}. 



(55) 
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Equation implies that 

H = R^HRT 1 ' 2 * = (R^HR-^V, (56) 

which means that if R 1 ' 2 were real, then H would be hermitian. However, this is generically 
not the case. A weaker result that will allow us to transpose easily most of the structure of 
the SUSY-QM case will be shown below: R is positive-definite - and hence R 1 ^ 2 real - when 
restricted to the eigenspace of the lowest eigenvectors. 



B. Gaussian development 

In the standard SUSY-QM case, it is useful to work in the basis in which H is Hermitian, 
mainly because in that basis the low-lying eigenvectors peak on saddle-points of the potential. 
It would seem that the analogous thing to do here is to go to an intermediate basis via 
R 1 / 2 . As R 1 / 2 is in general non- Hermitian, H is not Hermitian in this basis, only (complex) 
symmetric (H* = H^). We shall instead introduce a different basis \if) hR ) = e^\ip R ) and 
show that the \ip hR ) whose eigenvalues go to zero with the temperature are finite- variance 
Gaussians. Let us compute: 

fin 0H 

H = e—He~ — 

d 2 7 7 2 dV d d d d 2 V t t f ( 5? ) 

dp 2 2 4T % dpi dqt dqt dqidqj 

Although H' is not Hermitian, it will be more tractable than the original one. For each 
saddle point we propose a WKB form for the lowest eigenvectors: 

\TP hR ) = \rt R (q,p))®\^ h f R ) , . 

58 

In appendix O we show that the Gaussian so defined has finite variance. This is because 
the matrix B in ()58)1 is not singular, as it is for unstable saddle points in the original 
basis (reflecting the fact that in this basis the lowest eigenfunctions are not concentrated on 
saddles, see Section IVlII|) . We conclude that eigenvectors having eigenvalue zero (to leading 
order in T) are in this basis, just as in the ordinary SUSY-QM case, Gaussians peaked on 
saddle points: those with k fermions on saddles of index k. In this case, however, right and 
left eigenvectors do not coincide. 
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C. Zero-eigenvalue subspace 

In ordinary SUSY-QM, eigenvectors with exactly zero eigenvalue are annihilated by both 
Q and Q. This is evident in the basis in which H is Hermitian and Q the Hermitian conjugate 
of Q. In the present case, the proof is slightly more complicated. We start by writing: 

RH = Q ] RQ + Q ] RQ. (59) 

Clearly, the vectors annihilated by H and by RH are the same, and they must satisfy: 

{^\RH\^) = {i>\Q*RQ\^) + MQtRQty) = (Q^\R\Qi/j) + = 0. (60) 

If R were positive definite, then this would immediately imply that both Q\ip) and Q\ip) are 
zero. In fact, as mentioned above, R is not positive-definite, but one can show that it is so 
when restricted to the subspace of eigenvectors whose eigenvalues vanish to leading order in 
T (a space which obviously contains the eigenstates whose eigenvalues vanish exactly). To 
show this, we develop \ip) as in the previous subsection: 

|0> = ^a c e-^ c >®|r>, (61) 

c 

where the sum runs over the critical points 'c' of V(q). \g c ) is the normalized Gaussian 
centered on the critical point while \f c ) is the corresponding (normalized) fermionic part. 
Let us compute (<j)\R\<p): 

(<P\R\<P) = ^« c a c '*(e-^ c ®/ c |e^PJ|e-^/®r') = ^a c a c '*(/|P|/)(r| J\f c ). 

c,c' i,j 

(62) 

Two Gaussians centered on two different critical points do not overlap in the small temper- 
ature limit: 

(g c \P\g c ) = b c ^C c with C c > 0, (63) 
and the scalar product reduces to: 

(0| J Ri0) = ^i«fc c (riJir)- (64) 

c 

It remains to show that the fermion contribution is positive. We can assume that we have 
diagonalized the matrix of second derivatives of the potential at each saddle point, and for 
each direction i we have (see Appendix (0) that either Vu > 0, and then 

(fi\J\n) = 1, (65) 
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or Va < and 



(fiWi) = 2[7 2 (7 + vV - 4V«) - 2V«(27 + vV - 4V«)] > 0. (66) 

We have hence shown that in all cases within this subspace either 

(<j)\R\(j)} > or |0) = 0. (67) 

which, when applied to (jSOjl implies that eigenvectors with exactly zero eigenvalue are un- 
paired, just as in SUSY-QM. 

Paying the price of a loss of generality, this WKB approach shows in a more intuitive 
way the organization of the spectrum below the gap. The eigenstates can be seen - in an 
intermediate basis - as Gaussians centered on critical points. This Gaussian development 
enables us to show that there is no pairing in the zero eigenvalue eigenspace. 

VI. TIME DEPENDENT SUSY 

In this section we study the Kramers problem with a periodic time-dependent Hamil- 
tonian Ti,(q,p,t) of period r. We first introduce the Floquet formalism and extend the 
supersymmetry of the time-independent case to the time-dependent one. We then briefly 
show how this formalism allows us to rederive the Lefschetz formula and even prove the 
strong Morse inequalities when the friction is strong enough to prevent the proliferation of 
orbits of long periods. 

A. Generalized operators 

The Fokker-Planck equation associated with the Kramers system can be written: 

^(q,P\t)) = -H K (q,p;t)\ij(q,p;t)}. (68) 

Floquet theory is based on proposing solutions through the ansatz: 

\i>(q,p;t)) = \u(q,p;t))e- xt , (69) 

where \u(q,p; t)) is periodic of period r and the imaginary part of A can be chosen in the 
first 'Brillouin zone' [—-,-]. Equation (|68|) then becomes: 

(H K + ^) \u(q,p;t)) = \\u(q,p;t)). (70) 
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An alternative way to introduce the Floquet representation |17jj is to start from the 
stochastic equation. Introducing a variable 9 which grows linearly in time, we can write 
the Langevin equation (JHJ) as a system evolving with the time-independent Hamiltonian 
H(q,p,9): 



&H 

dpi 

% 
6 = 1, 



7— + ^qfrji 
dpi 



(71) 



with 



0(0) = 0. 

The equation (fTOj) leads us to a Kramers operator in the space (q,p, 9 

'dU d dU d 



N 



A 



N 



1=1 



d_ 

w 



(72) 



(73) 



-f dpi V dpi dpi J ~ \ % dpi dpi dqi 
To construct the supersymmetry of this operator, we first generalize the operator Q by 
introducing new fermion operators do and a\: 



Q - ia e 



d_ 
89 



N 

E 

i=l 











dqi 1 dpi 



bi 



d_ 
89 



(74) 



We can then extend Mr by adding fermion creation and annihilation operators: 
EI = M K + 5 5 6|Oj + 7 - oiftj + (jbhi + b]bj - a]di) 



+ 



dqidqj 

d 2 n 



dpidpj 



+ 7" 



d 2 n 



dpidpj ' 
d 2 H 



dqidpj 



(75) 



-aW 



d qi 39 'dp t d9j 1 d Pi d9 
HI commutes with Q, and we can also construct the other generator of supersymmetry: 



—i 







7 







i on 



i on. 



d 1 &H, 



dqi dpi T dqi T dp. 



— a 



satisfying: 

Q 2 = Q 2 = 0, T[ 
Let us also note that ae is a symmetry: 



dpi T dpi 



M. 



0. 



(76) 



(77) 



(78) 



implying that if is a right eigenvector of H, then ao\ip R ) is a degenerate one, or zero. 
In the following sections, we will use the supersymmetry to study the eigenvectors and use 
it to derive relations between the number of periodic trajectories. 
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B. Structure of the spectrum — Quadruplets 



We shall prove here that the right eigenvectors are in fact duplicated in the following 
way: one can build a basis {|^ )} = {\4>i), \Xi)> \Pi)} °f the space annihilated by ag which 
satisfy (|28|) and complete it by a degenerate free family {l^ 61 )} = {\<f4), \xt)i \pf)}i having 
non zero component along a\. This basis satisfies (see figure |3}: 

a e \ X ° l ) = a e \$) = a fl |p?) = 

a*lxf> ~ lx?> ~ 10°) a e \p^) ~ 

?> = lx?> QIp°) = o v|v)gw^k°> 

°> = \ti) = \xi) m) = \p") 





FIG. 3: Structure of the spectrum of EL The vertical axis represents the different complex eigen- 
values. The horizontal axes represent the number of fermions of type and the total number of 
fermions respectively. One can find a basis of the eigenvectors annihilated by ag which is organised 
as in the time independant case. It can be completed by a degenerate linearly independent family 
to build a basis of the whole space. 



To construct a basis as (|79|L one first remarks that the procedure followed in section Hi C 21 
is still valid in the subspace annihilated by ag. One can thus construct a basis {|p°), \4>T)i \Xi)} 
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satisfying (|25|l . Next, let us define {l^ 6 )} as: 

|tf) = Q|$) = Q|0°> + ioJ|^> 

|X?> = Q|0f> = (QQ + ^M) - ~aJ| X °) (80) 

Ip^QIp°> = QIp°> + ^4Ip°>- 

a# sends {|"0 )} back to {IV' )} ( U P to constant factors), which proves that {I'*/' 9 )} is an 
independent family. The whole family is obviously also independent, it can generate any 
vector annihilated by ae, thanks to the family {|^ )}, in particular it can generate the family 
(I?/; ), |a^ )} which is a basis of the whole space. 

Let us look at the (|p°)}- They are annihilated by Q but not the image by Q of any 
other eigenvector. The dimension of the space generated by such eigenvectors in the k 
fermion sector is the kth Betti number of the phase space {p, q} (see section Hi C 2j) . The 
dimension of the space generated by {|pf)} in the k fermion sector is then equal to Bk-i 
(see figure El). 

In the following, we will call "paired" states and "unpaired" states, the eigenvectors gen- 
erated by {\xi), \Xi), |0i>} and {|p°), \pf)}, respectively. 

As in any Floquet problem, the spectrum is organized in Brillouin zones. This can be seen 
directly as follows. Consider the family of periodic operators O m = e 2mme / T with m integer. 
Clearly, [H, O m ] = 27rimO m , which implies that if \ip R ) is an eigenvector with eigenvalue A 
then O m \ip R ) is an eigenvector with eigenvalue A + 2-irim. This redundancy is eliminated if 
we consider the operator e _nrH restricted to a particular starting point for 9: 

Tr ((9 = nr \e~ n ™\ 9 = 0» |° thcr = Tr [Te^M*] . (81) 

where 'other' means that the trace is taken over all variables (including fermions) except 9, 
and H{t) is: 

d / d dU &H\ d dU d 2 H L+ d 2 H _ t _ 

Hit) = -— 7 T— + — + 7 — + — — + -r^^blaj + 1 ——b\b 



dpi \ dpi ' dqi ' ' dp^J ' % dpi ' dq^qf 1 3 dpidpj 1 3 

d 2 n t , d 2 n . , t , t , t . / d 2 n d 2 n \ , t d 2 n t 

dp-dp- 4 ' 3 + dq-dp-^ + h & ~ ^ + [dq-dt + ^dp-dt ) blae - dp~dt aW 

(82) 
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C. Trace of the evolution operator 



Let us start by building a generating function just as in the previous section. We fix the 
starting point at 9 = and compute the trace: 



J^H{t)dt\ 



k fcrm 



(83) 



T(X,nr) = £A*Tr ((9 = nr | e - H | 9 = O))^ = £ ^ ^ 

k k 

where 'other' again means that the trace is over all the variables except 9. In the k fermion 
sector, it is divided in two parts: 



T k = Tr (V e~ ^ H ® dt ) = Tr ( T e" ^ H ® dt ) + Tr ( T e~ ^ m dt 

\ ' k fcrm. V / k fcrm. V 



paired 
k fcrm. 

(84) 

Let us show that unpaired eigenvectors have eigenvalue 0. Suppose that H\p° i ) = A|p°) with 
A not an integer multiple of i that is, A 7^ inside the first Brillouin zone. As Q and ag 
annihilate |p°): 

^b?> = QQb?> = ^). (85) 



Since Q — Q — idg-j^, can be written 



QQ\p°) 



(86) 



Because A is not zero inside the first Brillouin zone, O = A — J| acting on periodic functions 
of 9 is invertible. As [Q, O] = 0, one also has [Q, O" 1 } = 0. (j86|) then reads 



(87) 



which contradicts the fact that \p < f) is unpaired by Q. This shows that the unpaired eigen- 
states have eigenvalues zero, and the trace over them simply gives their number: 

unpaired 



Tr (je-irmdt 



k fcrm. 



B k + B k _i 



so that: 



2N 

T(X,nr) = (1 + A) ^ X k B k + ^ A fc Tr (r e~ ^ 



H(t) dt 



k=0 k 

On the other hand, because of the quartet structure, 



paired 
k fcrm. 



Tr (Te-K TH{t)dt ) = Tr [j e ~ ^ H{t) dt 

V J k fcrm. V 



paired without a' g 
k fcrm. 



-Tr (r e -W TH ® dt 



(88) 



(89) 



paired with 



k fcrm. 



(90) 
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leading to (see figure EJ): 

T(A,nr) = (1 + A) £ \ k B k + (1 + A) £ A fe Tr (V e" ^ H(t) * 

fc=0 k 
2N 2N 

T(A, nr) = (1 + A) £ \ k B k + (1 + A) £ A fc (i? fc (nr) + R k+l {nr)). 



without a# 
fe fcrm. 

27V 27V V ; 



fc=0 fc=0 

-nrH 



where again we have denoted Rk(nr) the partial trace of e~ n over fc-fermion states anni- 
hilated by ag and not by Q. 

We now perform the low temperature 'semi-classical' evaluation of the trace. Because we 
only need the trace restricted to states without ag, all the calculation in section III carries 
through without modifications, since the last two terms in (|82|) that are absent in (|12|) vanish 
in this subspace. We have then: 



T(A,nr) = ^A fe Tr (je^ H ® dt ) 

k 

= (l + A)^A fc Tr (VeJT^) 

k 



k fcrm 

without6> 



k ferm 

det(l + XU c (nr)) 



(92) 



v ' ^ |det(l — f/ c (nr))| 

noiseless v 
orbits c 

so that, finally: 

IN 2N 



E t'l^Ztl = E + E W + R^inr)). (93) 



Idetfl - U c (nr), t 

noiseless 1 v V //I fc=0 fc =0 

orbits c 

For A = — 1, this gives the Lefschetz formula. For any n: 

2N 

£ sign(det (1 - f/ c (nr))) = £)(-l)*B fc . (94) 

noiseless fc=0 
orbits c 

We can also obtain strong Morse inequalities for the case in which the total number of 
orbits is finite. This will happen if the friction is sufficiently strong and the time dependent 
forces are sufficiently weak. 

Consider first the case in which there are only orbits of period r (e.g. the case of a 
dissipative, adiabatic evolution). We concentrate on an orbit and its corresponding U c {r) 
with eigenvalues u^, Ui 2N . We order them so that u^.-.u^, verify |itj. | > 1 and Ui r+1 , Ui 2N 
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verify \ < 1, by definition the Morse index of the trajectory is then equal to r. We compute 
for large n the l.h.s. of (J9*H|) : 



det(l 4- WW) _ E k E „ y. , y, = 



|det(l -U c (nr))\ ^ ^ \1 - 

1 V V ,n k h,...,j k IU=1 I 1 "il fc Jl,.»Jfc 1 41 

since the only contribution which survives corresponds to {ji, ...,jk} = {ii, ■ ■, i r } 
In the limit n — ■> oo, (|93|) then becomes: 



2Af 

fc=0 fc fe 



A fc fi fc + A fe (^. + i(oo) + i? fe (oo)) = * kM k, (96) 



where is the number of orbits with Morse index k, and Rk(oo) the number of eigenstates 
of EI having zero eigenvalue (to this order in T) annihilated by ag and not by Q. 

Although we know of no concrete example, let us now outline how equation (|96J) would 
be derived for a system with a finite number of periodic orbits of several different periods. 
Consider again the limit of large n, but taken for n prime. The path-integral evaluation 
of the trace tells us that we have to sum the contributions of all orbits of period nr, that 
is, the repetitions of n times the orbit of period r. On the other hand, as we have seen 
above, the spectrum of Te~ ^ H i t ) dt contains, in the presence of orbits of prime period = pr, 
multiplets proportional (to this order) to a number times the p different roots of unity, and 
their contribution disappears from the trace over n cycles, since n is not a multiple of p. 
Hence, formula (J96|) is still valid, but now the Rk{oo) counts only the eigenstates of EI having 
zero eigenvalue (to this order in T) annihilated by ag and not by Q that are not part of a 
multiplet. 

The fact that the Rk(oo) are positive integers in equation (J96j) constitute the strong 
Morse inequalities, valid for relatively large 7 and/or small intensity of the time-dependent 
potential so that there is a finite total number of periodic orbits. 



VII. STOCHASTIC DYNAMICS IN THE fc-FERMION SECTOR 

In this section we describe a stochastic dynamics that corresponds to the extension of the 
Kramers equation to evolution of vectors with k fermions. We restrict ourselves to the case 
in which Ti is time- independent, although the generalization is straightforward. The purpose 
of this exercise is twofold: first, as we shall discuss below, we intend to use this dynamics as a 
practical method to find reaction paths and other structures in phase-space, and second, we 
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shall use the equations obtained to construct explicitly the low-temperature eigenvectors in 
the /c-fermion subspace, thus completing Morse Theory. In fact, the derivation is practically 
identical to the one for SUSY-QM in (J. 

We wish to find a stochastic process such that it somehow represents 

^ k {x,t)) = -H\il> k {x,t)), (97) 

where \ipk(x,t)) has k fermions. For large times, the \ipk(x,t)) will be a combination of 
states whose eigenvalues have small real parts. 

For zero fermions the dynamics associated with ()97|) is clearly Hamiltonian + noise + 
friction. Let us discuss the one-fermion sector in some detail. A one-fermion wavefunction 
has the form ipi(x)c\\—) , in phase-space variables. Equation (|97|) reads, for the components 

d 

— ipi(x, t) = -H K if)i{x, t) - Aijipjix, t). (98) 

Consider first one particle with a 2iV-component normalized vector u attached to it. The 
position of the particle evolves as a Langevin process and the vector u as: 

iii = —AijUj + N(u)ui, (99) 

where N(u) = J2ki u kAkiUi enforces the constancy of the norm Yli' 11 ! = 1- The only effect 
of the vector u on the dynamics is that we further impose that each particle has a creation- 
annihilation average rate = —N(u). From @ and (|99|h we have that the joint distribution 
function !F(x,u, t) evolves then as: 



~dt 



-H K - N{u) + — ( A H u j - N(u) Ul 



du 

One can check, using integration by parts, that 



T. (100) 



if>i(q,t) = J d\ Ui T(q,u,t) (101) 

will evolve according to 1)98)1 . 

The /c-fermion generalization is straightforward. The dynamics (197)1 for a vector 

^= E i> il ,..., ik (x)4 l ..4 k \-), (102) 

ii,...,tfc 
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where the V'u,...,ik( aj ) are totally antisymmetric, reads, in components: 

i>ii,...,i k = -H K Ai,...,i k ~ Y Y (-l) n(<T ' a) Ar(ii),a ^(fa),...,*,...,^). (103) 

where a denotes all permutations of k indices, and n(a, a) is the sign of the permutation 
(z'l, i 2 , ■ ■ ■ , a, . . . , ifc) — > (a, cr(ii), cr(z 2 ), • • • , cr(i fc )). Again, the particles follow equation (j9j), 
while the equations of motion for the v read: 

= - Y Y (-l) n(a ' Q) Ar(ii),a ««r(i 3 ),...,a > ...,a(i fc ) + A/"(t>), (104) 

with 

AT(t>) = Y v ^ Y Y(-v n{a ' a) A ^ ^ 2 ),..,<v..,^ fe )> ( 105 ) 

ii,...,ifc cr a 

thus preserving the normalization £\ i v^ ^. As before, there is cloning with rate J\f(v). 
It is easy to see that average of v^,...,^ indeed evolves as the V'ii ,.»>** 0*0 ■ 

The equations of motion for (v, J\f(v)) give the expansion rate of a volume element driven 
by the dynamics. Given a point x, we can write a volume element around it as: 

V k = 8x x A 5x 2 A ■ ■ • A 5x k = M Y v h,-,ik^h A ■ ■ • A ej fc , (106) 

where A is the external (wedge) product, are the basis vectors and the form fj lv ..,i fc is 
normalized. The orientation and norm of the volume element evolve along a trajectory of 
the particle. Equation ()104|) gives the evolution of the orientation v, and M — A4 gives the 
expansion rate. 

Low-lying eigenstates of H and geometric structures 

The right eigenvectors (x\ip^) with 'low' eigenvalues — with real part going to zero in 
the limit T — > — are concentrated on the following structures: 

• For k = 0: local minima. 

• For k — 1: paths originating in saddles of index one, spiralling down (thanks to 
friction) to local minima. 

• For any k: the /c-dimensional surface generated by all paths decreasing in energy 
emanating from a saddle of index k. 
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The argument is the same as in the purely dissipative SUSY-QM case Q: we consider that 
the surface descending from the saddle of index k is uniformly covered with particles with 
their fc-form v at each point tangential to such a surface. Both features are preserved by 
the evolution, as seen in the previous subsection. First, as the particles go downhill, the 
v attached to them change so as to remain tangential: this is because their evolution are 
precisely based on the linearized evolution on the tangent space. Secondly, the cloning 
rate matches exactly the expansion rate of a small volume advected downhill. Hence, the 
distribution of particles descending, and the average value of the forms fij,...,^ attached to 
them is left invariant by the T = dynamics. 



VIII. TRANSITION PATHS. 'REDUCED CURRENT' 



Part of the motivation for this work has been to use the higher fermion suspaces to find 
useful information on phase-space. In particular, in the Fokker-Planck SUSY-QM case, the 
low-lying eigenvectors of the one fermion subspace yield the transition currents between 
states. Here, as we shall see, the formalism itself tells us that, in order to compute the 
transition path, the relevant quantity is a 'reduced current', rather than the usual one. 
Consider first the Kramers equation (JTU1) : 

0P(q,p,t)_ D/ , T _ fdJ q . d,J n 



M - H " P ^ = - divJ --{-^ + ^)' < 107 > 

which defines the current: 

Inspired by the Langevin/Fokker-Planck case, we apply the operator Q to a distribution to 
obtain a current. We get: 

(-t)T QP(q,p,t) = \^ R ) = ,C4\-) + J T p fb\\-), (109) 

where we have defined the reduced current as: 

F cd = T , ^ dP(q,p) d _ ^ <9P(q, p) 

The reduced current has the following good properties: 
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• It differs from the current in a term without divergence, hence the fluxes over closed 
surfaces coincide. 

• It is zero in equilibrium. In a case with metastable states, it is small everywhere, while 
the true current is large within the states. 

It is easy to prove that Q is the only first-order differential operator that gives a 'current' 
field with the same fluxes as the standard current while being zero when applied to the 
equilibrium Gibbs state. 

In the Fokker-Planck SUSY-QM case, the supersymmetry allows us to obtain the usual 
transition current on the basis of 1-fermion low lying states. Here, the supersymmetric for- 
malism itself has suggested a new definition for the current, and a practical way to determine 
it (on the basis of simulating equation ([98)1 ). The reduced current is more relevant than the 
total one as it is directly related to passages rather than to phase space orbits within a state. 

An example. 

The example which we have studied is that of one particle evolving in the one- dimensional 
double-well potential V(q) = (q 2 — l) 2 . We have carried out simulations using the algorithm 
presented in the section IVHI and compared it with a determination of the current by direct 
simulation, see figures l4l and Io^3o1 | . 




FIG. 4: Reduced current for the double-well potential of section IVIIII One can see that the 
structure is concentrated on two damped paths, starting in the saddle point, and spiralling down 
to the two minima, respectively. The width of the structure is ~ \J^fT . 
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FIG. 5: Same as figure |IJ with smaller damping 
IX. CONCLUSION AND PERSPECTIVES. 

We have shown how the supersymmetry associated with the Kramers equation can be used 
to extract results of Morse theory in an elementary way. The strategy is very close to the one 
followed for the SUSY-quantum mechanics case, although somewhat complicated by the fact 
that the operators are non-Hermitian. Compared to the case of pure Hamiltonian dynamics, 
we have here a term of friction oc 7, and a term of noise oc \fyT which we can take to zero in 
several ways. Thanks to them, the Hilbert space on which the evolution operator acts is well 
defined, and the spectra of the operators are discrete. The wavefunctions associated with 
the eigenvalues with real part close to zero are concentrated on the structures associated 
with Morse Theory. 

The same methods can be used for a periodically time- dependent system, using a super- 
symmetric structure acting in the Floquet representation of the problem. The program for 
this case is however far from complete, below we mention some possible continuations. 

The advantage of this formulation is not only that it is entirely contained in (physics) 
undergraduate level, but that it makes a connection with situations of interest in physics 
and physical chemistry. Although we have used the low-temperature limit as a way to make 
the wavefunctions peak on saddles and other structures, other limiting situations yielding 
time scale-separation could have been invoked. For example, in many macroscopic systems, 
the thermodynamic states become mutually inaccessible (or almost), and thus play the same 
role as minima in the low-temperature case. The formalism applied here can then be used to 
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build a Morse Theory for states, their transition currents, and higher objects. The practical 
determination of metastable states and transition paths in this wider context is an active 
field of research, especially in complex energy landscapes relevant for chemical reactions. 
We have not dealt here with two elements that are needed in order to go beyond: 

• Degenerate Morse Theory, non-isolated orbits. 

• Proliferating orbits 

Degenerate solutions appear in a time-independent system as soon as we have to consider 
orbits rather than fixed points, corresponding to the freedom of choosing the starting point 
of the orbit. This degeneracy has already been discussed in the case of SUSY-QM 
and although it makes the treatment more cumbersome, it is a rather standard exercise in 
collective coordinates. Two cases when orbits arise naturally in a time-independent system 
are when forces do not derive from a global potential (like a magnetic monopole subjected 
to a magnetic field), or when the friction is scaled to zero with the temperature (e.g. 7/T 
finite). 

In systems for which the number of orbits grows exponentially with the period, a finer 
method of classifying orbits than the one used in this paper has to be put in place. This is 
the subject of Floer Theory, which we have only marginally touched. 

One last development that has been left out here is the case of Hamiltonian systems with 
no friction but with thermal noise (for such an infinite temperature situation, the phase- 
space has to be finite). Let us just remark here that in the corresponding low- noise limit, 
the supersymmetric formalism provides a method to study the separatrices and homoclinic 
orbits. 



APPENDIX A 



The generating function of the evolution operator can be written in path integral formal- 
ism as Lai: 



T(A, t) = Tr (\ F Te- 1 H ^') = £° V[q, p,n]5 (p i + 7 ^ + ^ 



s [ qi- 



rn 

dpi 



exp ("4^r / E 7 ^')^) 



(Al) 
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where W is defined by: 

W[q(t),p(t);t'} = f V[c,c}e~ Idts '( S ^^- lnX)+A ^ q] y\ (A2) 

and (,, ,) are & _ varices. U sing the F ou,e r _ tat ion of the , fenetien Q, 
one gets: 

T(A,t) = /%[fl 1 9 ) p ) p 1 ij]^ tf W^M^)-^ & ^^^]%p;i). 
•'go 

(A3) 

The question of factor ordering can be dealt with by choosing a convention in which the 
integral above represents the Kramers dynamics. In our case we need not worry about this, 
as we will only use the path integral as a bookkeeping device. The integration over the noise 
results in: 

T(A, t) = rP[ g ,g,p,p]e /d ^^(^^ + ^) + ^ 7T+ M^il)]H,(g,p ; t). (A4) 

Integration over p and q gives a 'Lagrangian' version of the path-integral: 

T(X,t) = J V[q,p}e~— Jdt ^^ + ^) s(q i --^-jW(q,p;t). (A5) 

The factor I/7T multiplying the action becomes large when the temperature goes to zero, 
and the path integral is then dominated by periodic orbits satisfying: 



dH OH 

dH 

dpi ' 



dqi dp 



(A6) 



One can linearize equation (|A4|) around such an orbit by putting X{ = x\ + VTx'^ Xi = 775?^ 
and consequently: 



m _ dH 

dxj dxj 



VTx', 



, d 2 n 



1 >p 



3 dxidxj 



(A7) 



One then gets for the contribution of the orbit to this order: 



T c (A,t) = / V[q\q,p\p}e 



J lq j a Pi a qj 



+Pj a qi e Pj 



q ,p q ,p 



Jdtqi (g<- 



8 Z H 



P J dpidpj H 3 dpidqj 

q ,p q 



^W(q c ,p c ;t) 



(A8) 

(the leading contribution to the action vanishes). Following our conventions for factor or- 
dering, one recognizes the path integral representation of the trace of the evolution operator 
associated with the time-dependent harmonic Hamiltonian ([31)1 . 
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APPENDIX B 

We now compute the trace corresponding to the bosonic degrees of freedom. Let us 
assume that the probability distribution is Gaussian: 



P(X,t) = exp f-1 - Xf(t)) (X j - X°(t)) + C{t)\ 



where Xf(t) has to be determined and C(t) is just a normalization factor. P(X,t) evolves 
with flUTJ): 



dP d 



dt dxk 



d 

Dki- h A c kj Xj ) P, 



fB2l 



which gives: 

ap 



(X, - - X°)(PPP - BA% - (Xi - X?)X]Bj& - DB kk + A% k P. (B3) 



dt 



On the other hand, differentiating directly (IB II) . we get: 



OP 



B 



-f(X t - Xf)(Xj - X°) + B y X°(X, - X°) + C{t) 
Equating (JB3J) and (IB4|) : 

-B = 2 (BDB) - (BA C ) - (BA C ) ] 
BX° = -BA C X° 
C(t) = —Tr(DB - A c ). 

The first equation implies for B '~ : 

—B~ l = 2D- AB 1 - B- l A\ 

dt 

which can be integrated to give: 

B 1 = 2 f\{t)U-\t')DU ] ~\t')U\t)dt' + U(t)a Q U\t). 
Jo 

Multiplying the second equation by P _1 on the left, one gets 

vO _ AC 

— ^Im^mi 

which means that X° follows the noiseless evolution. Equation (j3Tlj) tells us that 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 
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are solutions of (JB8J). 

From the normalization of P, we get: 



det B ® = e c(t) (B10) 

which satisfies the third equation of (jB5|) . Starting with P(X, 0) = 5(X — Y), i.e. X°(0) = 
Y and -By-(O) = lim A5y, the density is: 

A^oo 



P Y (X,t) = J ^^ e-^W-^WW-^y,). (Bll) 



S <3 



The trace now reads 



WOff^/^J ^ , (B12) 

' ( 2 ") N V det (1 -£/=(«))» det B det (1 - U<(t)) 

that is: 

/ P y (Y,t)dr = -—. (B13) 

J | det (1 — f/ c (t)) | 1 ; 

Therefore, the bosonic contribution of a classical periodic orbit is \ dct (i-u c )\ • 

The fermionic counterpart to this contribution is obtained by taking 

'EX'TrfTe-f^) = $> P E (-K-c^ e~^A^4 c \ |->, (B14) 

\ / p term, 

p p ii,...,t p 

where the time-order is defined along the noiseless trajectory. Let us first calculate the terms 
(-\c h ...Ci p T e _ ^ c<Cj ct / ..ct i |-). Using the fact that: 

(Te-J?^) c t (re/o*^) = . ( t ) 4, (B15) 

which can easily be seen by differentiating right and left with respect to t, and denoting 

= 7 e ~t A ii c * Cj , we get: 

(-\ Cll ...c tp o4 p ...cl i \-) = (-\c hl ...c ip o4o- 1 o4 p _ l o~ l ...o4 1 o- 1 \-) 

= l-\c- c U c ■ ct U c ■ ct ^U c -cU-) 



3v 3p-i n 



(B16) 



= En^H^,...^i,...ii-). 

jl~.jp k=l 

For a bracket to be non-zero, ji, ■ ■ ■ ,j p must be a permutation of i±, ■ • • , z p . The scalar 
product is then the sign of the permutation, and one has: 

(-\c tl ,...c lp 04 p ,...c\ i \-) = det p U^_ ip , (B17) 
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where det p £7? i is the minor of order p of U c (t) associated with the directions i±, ...,i p . 
We can now compute the generating function of the fermionic orbit's contribution: 

V A p Tr Te~ ti A « c * c > = V X p V det E/£ __ A 

— J p fermions J ' p ' ' p 

p p (B18) 

= det (1 + XU c (t)). 

Putting all together, the generating function of the contribution of a classical periodic 
orbit £ c' is: 

det (1 + XU c (t)) 

T c (X,t) = \ B19 

v ; det (1 - U c {t))\ v ' 



and we hence have: 



, . det(l + Af/ c t) 

T(X,t) = V \ j4rr. B20 

K J det (1 -U c (t)) V ; 

noiseless 1 \ v ' / 1 

orbits c 

APPENDIX C 

Let us consider a single saddle point (we drop the label 'c'), which we assume is at 
the origin. We develop V to second order in q i} and may assume further that (0) is 
diagonalized. We can hence treat each mode separately, so for ease of notation we drop the 
sub-indices We wish to construct for each mode of each saddle point, an eigenvector 
that is zero to this order: 

H c \^ h ) = 0, (CI) 

that is, for some A: 

(V'tfa + 7^6 - a ] b) \ij h f R ) = -X\^) R ). 

Let us first have a look at the fermionic part. In the basis |— ), \a*), \b<), \a*b*), the 
fermionic part of (jC2|) reads: 

(0 o o o\ 

o v" ~' o 1 ■ ( C3 ) 

7 J 

The spectrum and the corresponding eigenvectors are easily obtained: 

A = <-> \^f) = |-> 

A = - ± -^/^-AV" <-> \4> h f R ) = - 7 ± y/j2-4V»)\a!) + 2V"\tf) (C4) 
2 2 J 
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Regarding the bosonic part of ()C2|) . the Gaussian form of the eigenvector allows us to 
compute the l.h.s: 

(i 2 2 

(B pp - -) 7 + ([ - iBl - B pq f- + (V"B pq - ^if- 
24 1 1 (C5) 

+(-2jB pq B pp + V"B PP + B qq )^j \^ R ). 

The prefactor of the quadratic terms must be equal to zero. On the other hand, the Gaussian 
must be well normalized, i.e. the eigenvectors of B must be positive. After a tedious but 
straightforward calculation, two solutions are possible, depending on the sign of V": 

• if V" > 0, then B pp = |, B pq = and B qq = This corresponds to an eigenvalue 
for the bosonic part, and the corresponding fermionic part is thus the vacuum: 

l^>o) = e-^ (p2+yV) ®|-); (C6) 



if V" < 0, then B pp = ^ 1 \^ V \ B pq = ^ and B qq = -V" ^* AV " . This corresponds 

to an eigenvalue — § + — 2 ^ or bosonic part, which is compensated by the 

corresponding fermionic part (—7 — \f^ 2 — AV")\a)) + 2V"\ti i ): 



|^R <0 ) = e -*^^(p 2 -^V)-^>, ^ ^ ( _ 7 _ v / 7 2_ 4V , ) | a t ) + 2V r "|6 + >) • (C7) 
The same development for if t leads to: 

• if V" > 0, then 5 PP = |, S pg = and B qq = This corresponds to an eigenvalue 
for the bosonic part, and the corresponding fermionic part is thus the vacuum: 

l^»>o) = e-^ (p2+yV) ®|-); (C8) 



if V" < 0, then B pp = ^^ 2ry 4V " , B pq = and 5 OT = -V"' ^ 1 \^ V " ■ This corresponds 
to an eigenvalue — ^ + — ~ 4V , for the bosonic part, which is compensated by the 



corresponding fermionic part (7 + \f^ 2 — AV")\a)) + 2\tf): 



|^ ;<0 ) = e -^#^(f 2 -^"9 2 )+^r^ ^ 7 + ^/ 7 2_ 4 y//)| a t) + 2 |6t)) . (C9) 



The structure of the spectrum can be directly seen in the eigenvectors: a stable direction 
corresponds to the fermionic vacuum while an unstable direction corresponds to a one- 
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fermion state. We have also shown that both the right and left eigenvectors with zero 
eigenvalue are Gaussians to this order in the new basis. 



J. Milnor, Morse Theory (Princeton University Press, 1963). 
E. Witten, J. Diff. Geom 17, 661 (1982). 
H. Nicolai, Nucl. Phys. B 176, 419 (1980). 
H. Nicolai, Phys. Lett. B 89, 341 (1980). 

G. Parisi and N. Sourlas, Nucl. Phys. B206, 321 (1982). 
S. Cecotti and L. Girardello, Ann. Phys. (N.Y.) 145, 81 (1983). 

S. Tanase-Nicola and J. Kurchan, Journal of Statistical Physics 116, Nos. 5/6 (2004). 
S. Tanase-Nicola and J. Kurchan, Phys. Rev. Lett. 91, 188302 (2003). 
E. Gozzi, Progress of Theoretical Physics Supplement 111, 115 (1993). 
E. Deotto, E. Gozzi, and D. Mauro, J. Math. Phys 44, 5902 (2003). 
E. Deotto, E. Gozzi, and D. Mauro, J. Math. Phys 44, 5937 (2003). 
A. J. Niemi, Phys. Lett. B 355, 501 (1995). 

A. J. Niemi and P. Pasanen, Phys. Lett. B 386, 123 (1996). 
M. Miettinen and A. J. Niemi, Phys. Lett. B 461, 89 (1999). 

H. Risken, The Fokker-Plank equation (Springer, 1996). 
M. Nakahara, Geometry, Topology and Physics (Adam Hilber, 1990). 
P. Jung and P. Hanggi, Phys. Rev. A. 44, 8032 (1991). 

J. Zinn-Justin, Quantum Field Theory and Critical Phenomena (Oxford Science Publication, 
1996). 

O. Cepas and J. Kurchan, cond-mat/9706296 (1997). 
H. Kleinert and S. V. Shabanov, Physics Letters A 235, 105 (1997). 

B. Helffer and F. Nier, Hypoelliptic Estimates and Spectral Theory for Fokker-Planck Operators 
and Witten Laplacians, vol. 1862 of Lecture Notes in Mathematics (Springer, 2005). 
R. Mannella, Phys. Rev. E 69, 041107 (2004). 

for a discussion on the different possibilities to connect a classical evolution with a thermal 

n 
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n 

[24] These symmetries have already been obtained by Kleinert and co workers [211, for an action 
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corresponding to an evolution of the form JJ) in the (Lagrangian) path-integral formalism. 
[25] Diagonalizability is not entirely obvious even in the pure Kramers (zero-fermion) case, see \ 
for a discussion. 

[26] For a time-independent system with non-conservative forces for example: a single periodic 
(non-constant) trajectory is in fact a continuous family of trajectories, each one being the 
same orbit with a different starting point. 

[27] For instance if the points where & = are isolated in pj's direction, then this assertion holds. 
This is for example the case of TL = (p — f(q)) 2 + V(q). In more exotic cases, for instance if 
the phase space is periodic mp and TL = f(q) g(p) + h(q) with f(qo) = f'ilo) = 0, the solution 
corresponding to q = qo and p = h'(qo) t is a periodic orbit, which is not localized on a critical 
points. 

[28] In fact, this symmetry holds generally for systems with energy symmetric under the reflection 

of p: H(q, — p) = 7Y(q, p) and conservative forces. 
[29] We are assuming again that the number of real eigenvalues Uj < —1 is even, otherwise the 

contribution is — A r . 

[30] The numerical simulation of the equation (|T]l. especially in the low friction limit, should be 
implemented carefully, as a finite integration step will produce large errors. We have used 
quasi simplectic integrators (see |22| and references therein) which provide reliable results in 



the whole range of friction we studied. 



